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ABSTRACT 

We present a new theoretical estimate for the birthrate of R Coronae Borealis (RCB) stars that is 
in agreement with recent observational data. We find the current Galactic birthrate of RCB stars to 
be « 25% of the Galactic rate of Type la supernovae, assuming that RCB stars are formed through 
the merger of carbon-oxygen and helium-rich white dwarfs. Our new RCB birthrate (1.8 x 10“^ 
is a factor of 10 lower than previous theoretical estimates. This results in roughly 180-540 RCB stars 
in the Galaxy, depending on the RCB lifetime. From the theoretical and observational estimates, 
we calculate the total dust production from RCB stars and compare this rate to dust production 
from novae and born-again asymptotic giant branch (AGB) stars. We find that the amount of dust 
produced by RCB stars is comparable to the amounts produced by novae or born-again post-AGB 
stars, indicating that these merger objects are a viable source of carbonaceous pre-solar grains in 
the Galaxy. There are graphite grains with carbon and oxygen isotopic ratios consistent with the 
observed composition of RCB stars, adding weight to the suggestion that these rare objects are a 
source of stardust grains. 

Subject headings: ISM: abundances — novae, cataclysmic variables — binaries: close — stars: carbon 
— white dwarfs 


1. INTRODUCTION 


R Coronae Borealis (RCB) stars are a rare class of 
hydrogen-deficient supergiants with atmospheres com¬ 
posed primarily of helium and carbon, where the ratio 
of c arbon atoms t o oxygen atoms exceeds unity (C/0 
>1, Clayton 19961. They also show abrupt changes to 


their apparent visual magnitude. RCB stars are some¬ 
times observed at their maximum light but they have also 
been observed during their decline phase, in some cases 
showing sudden declines of up to 8 magnitudes in a few 
weeks. The star then typically stays faint for a period of 
several months ( Clayton|2012 1 before eventually recover¬ 
ing back to maximum brightness over a period of several 
months or longer. The steep decline in brightness is due 
to dust clouds which are ejected from the star and cover 
the photosphere when they are in the line of sight. Ob¬ 
servations indicate that this dust is made of amorphous 
carbon grain s which cover a range of sizes up to tenths 
of a micron (Garcia-Hernandez et al.||2011 2013 Jeffers 
et al.||2012D. 

The num ber of RGB stars kno wn in the Galaxy is 
about 100 (Tisserand et al. 2013), and more are being 
discovered or will be discovered in on-going and upcom¬ 
ing surveys of the Galaxy (e.g., WISE, GAIA, Euclid, 
LSST and SkyMapper). It is relatively easy to find RCB 
stars because they are variable and are bright, with ab¬ 
solute magnitudes between —5.2 < M y < —3.4, wher e 
the faint limit can be extended to —2.6 ( Tisserand|2012|) . 
To identify the object as RCB the typical decline in the 
light curve has to be identified but the process is sim- 
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plified because the minimum is very conspicuous in the 
light curve. Previous estimates suggest that the num¬ 
ber o f RCB stars in the Galaxy is higher, at over a 1000 
(e.g., |Glayton||2012|. In p|we discuss the observations, 
formation scenarios and theoretical models of RCB stars 
in more detail. 

Given that the known number of Galactic RGB stars 
has increased, we want to test whether these objects 
could be a viable source for producing carbonaceous dust 
in the Galaxy and if it would be possible to find their 
nucleosynthesis signature in pre-solar dust grains. Pre¬ 
solar grains are minerals that survived the formation of 
the Solar System and can be found in primitive mete¬ 
orites. Different stellar origins have been proposed for 
pre-solar grains including AGB stars, supernovae, no¬ 
vae, and post-AGB stars. The composition of pre-solar 
grains varies greatly and indicat es many different types 
of grains and stellar origin s (e.g.. 


Amari|[2005 Zinner|[2014 1. The bulk ot pre-solar grams 


Zinner 


1998 Lodders & 


both carbon and oxygen-rich types, are associated with 
the ejecta from supernovae or AGB stars. However, a 
few grains have be en associated with rarer sources such 
as nova outbursts (lA mari et al.||200i^ |Jos e et al.||2004|) 
and po st-AGB starsl Jadhav et al. 5013” Amari et al. 
(|2001b|) briefly discussed RGB stars as a source tor sil- 
icon carbide (SiC) grains of type A-l-B, which comprise 
about 4-5% of the measured SiC grains to date and are 
characterized b y their carbon isotope ratios, which show 
12 c/13c < 10. Amari et al. (2001b) discarded the idea 
because the majority ot carbon isotope measurements 
for RCB stars are ^^C/^^C > 40 as we discuss in 
Furthermore, SiC grains are not seen around RCB st ars 


(Lambert et al.||2001| Garci'a-Hernandez et al.||2011). 

In order to test these ideas we need to know the num- 
bers of RCB stars relative to novae and born-again post- 
AGB stars (see ®. However, of importance for dust 
production is theRust ejection mass and the lifetime so 
we also need estimates of these quantities, which we ob- 
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tain from theory or observations. This paper is organised 
as follows. We first provide an overview of the observa¬ 
tions and theory of RGB stars in ® In we use the 
results from binary population synthesis models to ob¬ 
tain the number of RGB stars expected from theoretical 
models where we assume that all RGB stars are formed 
by double white dwarf (WD) mergers. In 0we take ob¬ 
servational and theoretical estimates for the number of 
RGB stars to provide a dust-production rate. We also 
perform a similar calculation for novae and post-AGB 
stars. We finish in m with a discussion of our results. 

2. THE NATURE AND ORIGIN OF R COR BOR STARS 


There are two main competing theories for the forma¬ 
tion of RGB stars: 1) the final helium shell flash scenario 
or 2) the double degenerate WD merger scenario. The 
first scenario involves a WD central star of a planetary 
nebulae to undergo a final thermal pulse, which ignites 
the helium shell and cause s the star to expand to giant 
dimensions (Glayton 1996). While there is evidence for 
late thermal pulses in the general population of post- 


Duerbeck & Benetti 

1996[ [Glayton & De Marco[ 

TMtI 

Asplund et al. 

1997 

1, there are serious problems with 


this scenario expl aining the lifeti me of RGB stars. We 
refer to Table 3 in Glayton (2012) for an overview of the 
double degenerate versus hiial hash formation scenarios. 

The second scenario states that RGB stars are the 
merger product s of carbon-oxyg en (GO) WD and helium- 
rich WD pairs (Webbink 1984). The progenitor binary 
system experienced at least one phase of common enve¬ 
lope (GE) evolution (see Due to radiation of gravi¬ 
tational waves the system undergoes orbital decay which 
causes the two stars to approach each other and finally 
coalesce. During this merger, the less massive helium 
WD gets completely disrupted. The matter from the 
helium WD is believed to form a Keplerian disc which 
gets assimilated o nto the surface of t he accretor, where 
it starts to burn (|Jeffery et al. 2011). The remnant of 
the helium WD niight also f orm an extended en velope 
around this merging product ( Glayton et al.|[2MT ). 

RGB stars have pulsation periods on the order of 40 to 
100 days, from whic h pulsation mas ses of « 0.8 — 0.9Mq 
have been derived (Glayton 1996). Models have also 
shown that the merger of a GU WD and a helium WD 
is predicte d to result in a mass of about 0.96 ± O.ISMq 
( H an|1998 ), which is consistent with the pulsation masses 
( Glayton[I996 Glayton et aLpOlI ). It has also been hy- 
pothesized that such mergers may lead to thermonuclear 
explosions such as Type la supernovae ( SNe la), but this 
remains a topi c of in vestigation ( see e.g. Dan et al.|2014 


Pakmor et al.|[20T3 ). Han (1998) theoretically predicted 
the birthrate of RGB stars to be 1.8 x 10“^ yr“^. In 
w e update the predict ed RGB birthrate for the Galaxy. 

I Jeffery et al. (|2011|) presented surface abundances of 
RGB with respect to evolutionary models which as¬ 
sume that RGB stars are formed by a double degener¬ 
ate merger of a GO WD and a helium WD. The mod¬ 


els were compared to the observations (Asplund et al. 


2000 

whic. 


Rao & Lambert 2003 2008 |Pande^ et al. 2008), 


1 reveal that RGB stars are enhanced m several el- 


ements such as nitrogen, fluorine, sodium, aluminium, 
phosphorus, silicon, sulfur, and in some neutron capture 
elements produced by the slow neutron capture process. 


Observations also show unexpectedly low isotopic ratios 
of ^^O/^^O with ratios close to or less than unity ( Glay- 


ton et al.|2005 2007 Garcfa-Hernandez et al.|20i0 ). 'I'he 

ratios of ^^ G/^^G are generally high, with a lower limi t 
of 30 to 40 (Cottrell & Lambert|1982 Hema et al.|2012 ). 
There are exceptions, inclu ding V CrA which has a low 
12 c/13c ratio of » 3 to 4 ( R ao fc Lambert 2008), and 
VZ Sgr with 3-6 ( Hema et al.||20I2 r 


The elemental and isotopic abundances of RGB stars 
are reaso nably well matched b y merger models, such as 
thos e by Jeffery e t al. (|2011|), Longland et al. ( 2011]) 
and liV lenon et al.| (|2013|) . The models by|Jettery et al. 
(20111 assume most of tire observed surface composition 


of rGB stars comes from th e nucleosynthesis fro m the 
AGB phase of the CO WD. Jeffery et al. (2011) iden¬ 
tified the need for additional nucleosynthesis during the 
merger to explain the composition of some elements, such 
as fluo rine, which shows a bundances higher than pre¬ 
dicted. Menon et al. (2013) reproduces the abundances 
of C, In , U and 1' using a one-dimensional post-merger 
evolution and nucleosynthesis model based on realistic 
hydrodynamic merger progenitor models. However, the 
large Si and S abundances observed in the low metallicity 
RGB stars are not reproduced, indicating that some high 
temperature burning at T « 10®K may be occurring. 


3. RATES OF R COR BOR STARS FROM POPULATION 
SYNTHESIS 


Theoretical (or extrapolated) estimates for the number 


Lawson et al. 

1990 

Han 

00 

mate the num 

ber 0 

RGL 

stars 


2U12p . 'lb esti- 


RGB stars expected to currently be 
in the Galaxy, we calculate theoretical merger rates be¬ 
tween helium-rich and CO WDs (hereafter He-CO merg¬ 
ers) using the binary evolution population s ynthesis code 
StarTrack dBelczynski et al.(j2()02l |2008|). 

Recently, StarT'rack input physics have been up¬ 
dated to include a new prescription for accretion of 
helium-rich material on CO and ONe W Ds when sta¬ 
ble Roche lobe overflow is encountered (iRuiter et al.l 


2014). The updated prescription does not aflect double 
WDs that merge upon contact (like RGB progenitors), 
but does affect binaries involving He-CO systems that do 
not merge upon contact and thus lead to AM CVn sys¬ 
tems, ‘classical’ double-detonation Type la supernovae, 
or other transient e vents. We use the s imulations from 
the P-MDS model of Ruiter et al. (2014) to calculate the 
birthrate of RGB stars arising froiri He-CO mergers. All 
stars are evolved from the ZAMS, and the evolutionary 
histories are record ed until t = 13.7 Gyr is reached (see 
[Ruiter et al.| 20091 Section 2, for a description of initial 
separations, initial mass function, etc.). For these simu¬ 


lations we have adopted a ZAMS binary fraction of 70% 
and a burst star formation history, where the burst oc¬ 
curs at t = 0. A ‘burst’ simulation is the most powerful 
since it allows us to easily convolve our event rates (delay 
times; see below) with any star formation rate of choice 
to attain a realistic model star formation history. 

We note that the delay time distributiorj^ shape of He- 


® The delay time distribution (or DTD) is the distribution of 
events in time — this case, mergers - that occur following a hy¬ 
pothetical burst of star formation at t = 0. The DTD shape is 
a useful tool in quantifying progenitor ages for explosive stellar 
phenomena, e.g. Type la supernovae. 
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Figure 1. Total mass of RGB systems at time of merger grouped 
into mass bins of 0.05 Mq. We assume all RGB stars are formed 
through the merger of a GO WD and a helium-rich WD. 


CO mergers relative to CO-CO mergers is qualitatively 
similar starting from 500 Myr post-starburst, owing to 
the fact that similar physical processes play a role in de¬ 
termining the post-nuclear burning orbital evolution of 
double WDs. CO-He mergers do not occur with delay 
times < 500 Myr due to the low-mass nature of the sec¬ 
ondary star: A helium-WD progenitor takes longer to 
complete nuclear evolution than a CO-WD progenitor. 

Given the similar DTD shapes beyond 500 Myr, it is 
reasonable to assume a constant birthrate of RCB to CO¬ 
CO mergers at the current epoch (~11 Gyr). We find 
the number of He-CO mergers averaged over a Hubble 
time to be 2.65 x 10“^^ M 0 “^yr“^, where the mass rep¬ 
resents mass born in stars. This value is 25 % of our 
CO-CO merger rate: 1.06 x 10“^^ MQ~^yT~^. This rate 
is remarkably similar to the rate of Type la supernovae 
in Sbc-like galaxies such as the Milky Way: 1.1 x 10~^^ 
Mf 7 ^~^yr~^, or 7 x 10“^ SNe la yr“^ as cited by Badenes 
fc Maoz (2012) when a dopting a Milky Way stellar mass 
of 6.4 ± 0.6 X 10^° Mq ( McMillan||2011 ® We adopt this 
Galactic SN la rate in order to extrapolate our RGB pop¬ 
ulation to Galactic numbers: we take the specific rate of 
He-GO mergers to the observed specihc Galactic SN la 
rate which yields an RCB birthrate of 1.8 x 10“^ RCB 
stars yr“^. 

In Figure we show the distribution of total masses 
for He-CO mergers over a Hubble time calculated with 
StarTrack. The mass distribution peaks at « O.QMq, 
which is consistent with the theoretical estimates of lHanI 
(1998) as well as t he inferred RCB masses derived from 
pulsation studies (Clayton 1996 Clayton et al. 2011|). 
We note that in principle observationally-measured RCB 
masses are likely to be equal to or less than these theoret¬ 
ically predicted RCB ‘birth’ masses, because the merger 
product will likely lose some mass over time. 


® Note that our simulations are normalized p er unit mass born 
in sta rs where as the observational value from IBadenes &: Ma^ 
||2012|l is normalized per unit stellar mass in a given galaxy. We 
have not attempted to correct for this difference since it unneces¬ 
sarily introduces uncertainties into our rate. However, changing 
the normalization technique to more closely match that of obser¬ 
vations would lead to a slight increase in our rate, thus we expect 
even better agreement. 


We find that RCB stars are formed via three main 
evolutionary channels, all of which involve one CE event. 
Factors such as ZAMS mass and initial orbital separation 
play a role in how soon after star formation the He-CO 
merger occurs, as well as what the component masses and 
compositions will be at time of merger. Two of the evolu¬ 
tionary channels (hereafter channel 1 and channel 2 ) in¬ 
volve the merger of a CO WD and a ‘hybrid’ helium-rich 
WD, where the WD has a CO-rich core and a helium-rich 
mantle. Such HeCO hybrid WDs can be formed in cases 
where a red giant star is stripped of its envelope during 
binary interacti ons, and the helium core only undergoes 
partial burning (Tutukov & Yungelson||1996 ). The third 
evolutionary channel (channel 3j involves the merger of a 
CO WD and a He WD, the latter which never underwent 
any helium core burning. 

The main difference between channels 1 and 2 is that 
the ZAMS masses are typically larger by for both 

components in channel 1. Channel 1 CE events occur 
when the secondary is a red giant, and the WD merg¬ 
ers occur relatively quickly after star formation (within 
0.5—1 Gyr). For channel 2, the GE events occur while 
the secondary is an early AGB star, and the WD merg¬ 
ers occur 1—4 Gyr post-star formation. Ghannel 3 binary 
components are slightly less massive both on the ZAMS 
and at time of merger. The secondary star is found on 
the early AGB during the GE phase, and the WD merg¬ 
ers occur > 4 Gyr after star formation. 

We find the most common formation scenario to be of 
channel 2 type, followed closely by channel 3. In Eigure[^ 
we show a typical evolution for an RGB star formed via 
channel 2 from our population synthesis model. Other 
than birth on the ZAMS and the time of merger (marked 
i and x in the figure, respectively), we have depicted 
stages where either a stellar component first evolves off 
the main sequence (ii, v), initiates a mass transfer phase 
(hi, vi-viii), or becomes a white dwarf (iv, ix). Over 
the course of evolution, the binary undergoes a stable 
mass transfer episode when the primary star is on the 
Hertzsprung Gap, which follows into the red giant phase. 
After the primary star evolves into a He-rich (hybrid) 
WD, a GE event occurs while the secondary is an early 
AGB star. The CE leaves behind a naked helium-burning 
star in a relatively close orbit (1 Rq) with the He-rich 
WD. After the second star becomes a WD, it takes ~800 
Myr for the stars to be brought into contact under the 
influence of gravitational wave radiation. 


4. DUST PRODUCTION RATES 

An important factor when determining the dust pro¬ 
duction rates from various sources is not just the num¬ 
ber of them at any one time in the Galaxy, but also the 
amount of dust that each source produces. In this sec¬ 
tion we present estimates for the number of RCB stars, 
novae and born-again post-AGB stars and their dust- 
production rates. Table [l] is a summary of our results, 
where we highlight the source, the birthrate, the dust 
mass formed per event, and the dust production rate. 


4.1. R Coronae Borealis stars 

In the last decade the observe d number of RCB stars in 
the Galaxy has risen from ~40J Zaniewski et al.| 2005 ) to 
76 (jTisserand et al.||2013|). |Lawson et ai.| p99l)| extrap¬ 
olated from the nuniber of known RGB stars at the time 
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Table 1 

Birthrates, expelled dust mass and dust production from RGB stars, novae and born-again post-AGB stars in the Galaxy. 


Source 

Birthrate [yr 

Dust expelled per event [Mq] 

Dust production rate [Mq yr 

R Goronae Borealis stars 

1.8 x 10“® 

0.1 - 1 x 10-1 

1 - 50 x 10-1 

Novae explosions 

1 - 10 x IQi 

10-1° to 10-° 

CO 

1 

o 

1 

00 

1 

O 

Born-again post-AGB stars 

0.63 - 28 x 10-1 

4 x 10-1 

7.5 - 34 x 10-° 


t=o 

Ml=2.26 



a = 17.3 
M2=1.58 


t=813.4 

Ml=2.26 


o 


t=824.8 

Ml=2.25 



a = 16.6 
M2=1.58 


t=836.9 
Ml = 0.34 


o 

He-rich WD 


a=126.2 

M2=2.53 


t=1295.6 
Ml = 0.34 


o 



t=l468.4 
Ml = 0.34 


t=l468.4 
Ml = 0.34 



t=1469.9 
Ml = 0.34 


a=1.0 
M2=0.59 

CO WD 



t=2276.7 
Ml = 0.34 


o# 

total mass = 0.93 


MERGER 
M2 = 0.59 


Figure 2. Evolution showing one of the three principle evolution¬ 
ary channels through which RGB stars are expected to form, as 
predicted by our StarTrack model. Only the ten most interest¬ 
ing phases of evolution are shown (see text). Times are in Myr since 
birth on ZAMS, separations (a) are in Rq, and masses (Ml, M2) 
are in Mq. The phase of stellar evolution prior to WD formation 
is denoted as follows: MS = Main Sequence, HG = Hertzsprung 
Gap, AGB = Asymptotic Giant Branch, He = (stripped envelope) 
helium-burning star. 


to estimate that there are 200 to 1000 RGB (and cool 
hydrogen dehc ient carbon) stars in the Galaxy, whereas 


Clayton (2012) extrapolated from the number of known 
RGB stars in the LMC to arrive at a considerably higher 

number of 5,700 RGB stars. _ 

We perfo rm a similar extrapolation here. ITisseraiid] 


et al. p009 ) find 22 RGB stars in the Magellanic Clouds 
with l8 in the Large Magellanic Cloud (LMC) and 4 
in the Small Magellanic Cloud (SMC). The Milky Way 
Galaxy is roughly 10 times more massive than the LMC 
and 100 times more than the SMC. If we crudely assume 
that the number of objects scale with galaxy mass (P. R. 
Wood, private communication), we estimate that there 
are between 180 and 400 RGB stars in our Galaxy. Only 
part of the LMG has been searched by MACHO and 
OGLE and that the number of RGB stars in that galaxy 
may be significantly higher than 18. Either way, our new 
numbe rs are signihcantly less than estimated by |Clayton| 
(20121. Our new theoretical birthrate of 1.8 x 10“'^ RGB 
stars y r~^, a factor of 10 lower than predicted by Han 
(1998), results in « 180 — 540 Galactic sources when 
taking RGB lifetimes between 1 — 3 x 10^ years. This 
theoretical estimate is in good agreement with our ex¬ 
trapolated number of RGB stars. 

Observational estimates of the amount of carbon-rich 
dust elected from RGB stars ranges from 10~'^M^^ vr~ ^ 
(IClayton et al. 1992) to 1Q~^Mq yr“^ (Feast 1986). 
These numbers are consistent with the l atest observa¬ 
tional esti mate of 9 x 10~^Mf7^ yr“^ fromJClayton et al. 
(2011) and Jeffers et al. (2012). In Tablewe calculate 
our dust-production rates using \Q~^Mq yr“^; choosing 
the lower limit of Mq yr“^ would lower our dust 
production rates by an order of magnitude. Further¬ 
more, we have assumed that all RGB stars produce the 
same amount of carbon-rich dust, which may not always 
be the case. 

Taking the observed number of RGB stars in the 
Gala?^ to be ^ 100 we obtain a dust production rate 
of « lO^'^MQyr”^ , assuming a dust mass-loss 

rate of 10~^Mq yr“^. Using our newly calculated RGB 
birthrate, adopting the upper limit for the RGB lifetime 
of 3 X 10® years, and again assuming a dust mass-loss rate 
of 10“®Mq yr“^, our best theoretical estimate yields a 
dust production rate of « 5.4 x 10“"^ Mq yr“L 

In order to find a ha rd upper li mit, we take the theo¬ 
retical birthrate from Han (|l998 ) of 1.8 x 10“^ yr“^, an 


RGB lifetime of 3 x 10® years which results in 5400 RGB 
stars. This yields a dust production rate of 5 x 10 “®Mq 
yr“^. Our range of dust-production rates for RGB stars 
is summarized in Tabled 

One substantial uncertainty in this estimate is the life¬ 
time of RGB stars. That there are only about 100 found 
out of the « 200 — 5000 theoretically predicted may sim¬ 
ply reflect a sho rter lifetime for the RGB phase than es¬ 
timated by e.g., Glayton (2012). 
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4.2. Novae 

We now compare these rates to the dust production 
expected from novae. Similar to the situation for RGB 
stars, the frequencies and mass-loss rates for novae cover 
a wide range. Galactic nova rate estimates ha ve yielded 


consistent with estir nates of planetary n e bula birthrates 
of about 0.5 yr~^ by Zi.ilstra & Pottasch (1991), whereas 
Moe & De Marco (2006) hnd a higher value of 1.7 ± 0.3 


yr ^ for the entire Milky Way. Moe & De Marco (2006 1 
also give the post-AGB birthrate of 2.2 ±0.5 yr“^, which 
is estimated to be 90% of the white dwarf birthrate. The 


numbers from 10 yr ^ (Ciardullo et al. Ij90) to « 100 


yr ^ (Liller & Mayer 1987| Shatter 1997). Most esti- 


post-AGB birthrate from Moe & De Marco (2006) is 

- - 


mates are somewhere in the middle of this range at val- cl oser to our uppe r estimate of 3 yr 


ues of around « 30 — 40 yr 

^ (Hatano et al.|1997 

Nelson 

et al. 

2004 

Darnley et al. 

2006|). Nova rates are bet- 


ter constrained in M 31 than m the Milk y Way, with an 


estimate of 65 yr ^ ( Darnley et al.|[^006 ). 


The error estimates are highly uncertain. We take the 
upper limit of 100 yr“^ to ensure that we are not un¬ 
derestimating dust production from novae. Typical dust 
formation masses span a wide range from 10 “to 
an upper limit of « 1 O“^M 0 , where most of the dust 


produced by novae is carbon rich (Gehrz et al. 
We will take the lower and upper 


1998 
imits (see I'able 


I 


which yield nova dust production rates of M7°yf 
10 “® to IO-^Mq yr-i. 

Note that nova explosions do not always result in de- 
tectable dus t prod uction, with 3 out of the 25 sources in 
Gehrz et al. (1998) showing no dust, and another 4 with 
no information on the types of dust formed. 

4.3. Born-again post-AGB stars 

There is no expected dust production associated with 
the post-A GB phase because the stars are too hot to 
form dust (van Winckel 2003). The dust around post- 
AGB stars IS recycled AGB dust. The late and very late 
thermal pulse post-AGB stars are different: These ob¬ 
jects are born-again giant stars such as Sak urai’s Object 
that may pro duce their own, G-rich dust (Duerbeck 
BenettiHl996|). Born-agai n stars have been consid ered a 
source of pre-solar grains ( Jadhav et al.|2008 2013| ). Few 
born-again stars have been studied so we do not nave a 
good estimate of the observed mass-loss rates . Saku rai’s 
Object is the best studied and van Hoof et al. (2007) cal¬ 
culate a lower limit to the mass of the total ejecta to be 
6 X 1 O“^M 0 . The most recent estimate gives a dust mass 
of 4 X 1 O“®M 0 (A. Zijlstra, private communication). We 
assume that each born-again AGB star ejects the same 
amount of dust as Sakurai’s Object over the born-again 
lifetime. The lifetime of a born-again AGB stars is ap¬ 
proximately that of a thermal pulse {ttp ~ 10 ^ years), 
giving an average dust mass-loss rate of 4x 1 O“®M 0 yr“^. 

While we do not know the number of born-again AGB 
stars in the Galaxy, we can estimate ho w many there 
are usin g th e number of post-AG B stars. |Kamath et al.| 
(2014b) and Kamathet al. (2014a) recently estimated the 
number of post-AGB stars in the LMG and SMG to be 75 
and 34, respectively. These numbers include both single 
and binary post-AGB stars, and they derive a combined 
post-AGB lifetime of « 5 x 10^ y ears, consistent with 
lifetim es from evolutionary tracks (Vassiliadis & Wood 
1994). If we scale the number of post-AGB stars in the 
Magellanic Glouds to the Galaxy like we did for RGB 
stars we obtain between 750-3400 post-AGB stars. 

To calculate the birthrate of post-AGB stars in the 
Galaxy we require a lifetime, whic h has been theoreti¬ 
cally estimated to be « 10^ — 10"^ (Vassiliadis & Wood 
1994). We obtain 0.18 to 3.4 yr“^. The lower values are 


Iben (1984) and Renzini (1982) estimate that 10-25% 
of all stars that leave the AgB experience a late thermal 
pulse. Using our new extrapolation for the number of 
post-AGB stars in the Galaxy, and assuming that 25% 
of all post-AGB stars experience a late thermal pulse, 
we obtain between 188 to 850 born-again AGB stars. 
Using the expelled mass provided in Table we estimate 
that born-AGB stars produce « 7.5 — 34 x 

10-^Mq yr-i. 

Assuming a lifetime of 10^ years for a late thermal 
pulse allows us to estimate the birthrate of born-again 
AGB stars to be between 1.9-8.5 yr“^. The upper limit 
seems unreasonably high. Either we have overestimated 
the number of post-AGB stars or we have underestimated 
the lifetime of a late helium-shell flash. The low-mass 


(< 2Mq) AGB models of Karakas (2014) spend a few 
times 10^ years with a convective He-shell. In Table 
we assume the timescale of a He-shell flash is 300 years, 
which gives birthrates of 0.63 — 2.83 born-again post- 
AGB stars per year. These birthrates are much more 
reaso nable when compared t o the post-AGB birthrate 
from|Moe & De Marco (2006). 

Unlike RGB stars, not all of the dust from born-again 
AGB stars may be carbon rich but this is yet to be de¬ 
termined. Reducing the number of post-AGB stars that 
experience late thermal pulses to 10 % lowers the number 
of stars to 75-340, and the amount of dust produced to 
3-14 X 10 ^Mq yr Note that increasing the dura¬ 
tion of thermal pulse to 300 years will decrease the dust 
mass-loss rate. 

5. DISCUSSION AND CONCLUDING REMARKS 

We present a new theoretical birthrate for RGB stars of 
1.8 X 10“^ yr~^. a fac tor of 10 lower than the previous es- 


tima t e by Han (1998). We note that in the simulations of 
Han| (|1998|, it was assumed that all double degenerates 
that reach contact will eventually merge (Z. Han, pri¬ 
vate communication, 2015). Therefore, one would nat¬ 
urally expect a higher RGB birthrate than we obtain 
here, since in StarTrack we do not assume all double 
degenerates merge. Whether a double degenerate binary 
merges depends on the adopt ed mass transfer sta bility 
criteria for a given model (see Toonen et al. 2014) , and 
will st rongly impact the resulting merger rates. |Shen| 
(|2015|) recently suggested that every interacting double 
white dwarf binary may merge, which would increase the 
predicted number of RGB stars. Using our new birthrate 
we predict that there should be roughly 180-540 RGB 
stars in the Galaxy. We also estimate the number of 
RGB stars independently from the number known in the 
Magellanic Clouds, and arrive at an estimate consistent 
with our new theoretical predictions. 

We calculate the dust production rate from RGB stars, 
novae and born-again post-AGB stars with our results 
presented in Table We show that the dust mass-loss 
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rate from RGB stars is comparable to or higher than no¬ 
vae and born-again post-AGB stars, even when account¬ 
ing for the large uncertainty in the numbers of these ob¬ 
jects in the Galaxy today. Our simple arguments along 
with updated observational and theoretical estimates for 
the numbers of RGB stars in the Galaxy show that RGB 
stars are a viable source of carbonaceous pre-solar grains. 

Are there any pre-solar grains that show the chemical 
signature of RGB stars? Recall that RGB stars have high 
G/0 ratios, ^^G/ ^J^G ratios > 40 and ratios be- 

tween « 0.2 — 20 (Glayton et al.|20'07 Garcia-Hernandez 


et al.| 2010), and some show the sig nature ot hea.yy el¬ 


ement production by the s process (Jeffery et al. 2011 


Menon et al.||2013). Garcia-Hernandez et al. 12010) also 


attempt to derive ratios but only provide lower 

limits due to the non-detection of th e ^^G^^N line. _ 

The pre-solar grain database from [Hynes fc Gyngard] 
(20091 has a comprehensive list of all the measured iso- 
topic and elemental compositions of pre-solar grains from 
the literature. Here we focus on the graphite grains 
as these likely formed in a highly ^^G-rich environment. 
Furthermore, silicon carbide features have not been de- 
tected aroun d any RGB stars to date ([Garcia-Hernandez 


et al. 2011[ [2013 ). Amari et al. ( l993p observe that 


graphite grains with the largest excesses also have 
the highest ^^G/^^G ratios. These could be associated 
with an RGB origin. There are roughly 1800 graphite 
grains with carbon isotopic ratios that match RGB stars 
but only about 10 grains with measured ratios 

< 25. All 10 grains have large ^^G/^^G ratios (> 160) 
and high ^®A1/^^A1 ratios > 0.028, consistent with hy¬ 
drogen burning. |Menon et al. (2013) predict that Al can 
be produced, depending on tne interplay between the 
temperature and the burning timescale but they do not 
provide predictions for ^®A1 in particular. They also note 
that both Si and Ga is observed enhanced in RGB stars 
but their models do not reproduce these abun dances. 
Neutron cap tures also occur in the models by Menon| 


et al. (2013), which could explain the excesses of 
observed in the graphite grains. Unfortunately they do 
not provide isotopic predictions for comparison to the 
grain data. Future work is needed to further explore the 
viability of whether the grains with excesses in could 
indeed have originated from RGB stars. 
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